FROM CONSTITUENT QUARK TO HADRON STRUCTURE IN THE 
NEXT-to-LEADING ORDER: NUCLEON AND PION 
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We calculate the partonic structure of constituent quark in the Next-to-Leading Order for the 
first time. The structure of any hadron can be obtained thereafter using a convolution method. Such 
a procedure is used to generate the structure function of proton and pion in NLO. It is shown that 
while the constituent quark structure is generated purely perturbatively and accounts to most part of 
the hadronic structure, there is a few percent contributions coming from the nonperturbative sector 
in the hadronic structure. This contribution plays the key role in explaining the SU(2) symmetry 
breaking of the nucleon sea and the observed violation of Gottfried sum rule. These effects are 
calculated. Excellent agreement with data in a wide range of x = [10 -6 , 1] and Q 2 = [0.5, 5000] 
, GeV 2 is reached for proton structure function. We have also calculated Pion structure and compared 

it with the existing data. Again nice agreement is achieved. PACS Numbers: 13.60 Hb, 12.39.- 
x, 13.88 +e, 12.20.Fv 

<\ 

OO I I. INTRODUCTION 

^vq . Our knowledge of hadronic structure is based on the hadronic spectroscopy and the Deep Inelastic Scattering (DIS) 
' data. In the former picture quarks are massive particles and their bound states describe the static properties of 
f — , hadrons; while the interpretation of DIS data relies upon the quarks of QCD Lagrangian with very small mass. The 
hadronic structure in this picture is intimately connected with the presence of a large number of partons (quarks and 
gluons). The two types of quarks not only differ in mass, but also on other important properties, for example, the 
color charge of quark field in QCD Lagrangian is ill-defined and is not gauge invariant, reflecting the color of gluons in 
an interacting theory. On the other hand, color associated with a Constituent Quark (CQ) is a well defined entity. It 
has been recently shown that |lj Q one can perturbatively dress a QCD Lagrangian field to all orders and construct 
a CQ in conformity with the color confinement. From this point of view a CQ is defined as a quasi-particle emerging 
from the dressing of valence quark with gluons and q — q pairs in QCD. 



Of course, the concept of CQ as an intermediate step between the quarks of QCD Lagrangian and hadrons is not 
On! new. In fact, in the context of SU(6)XO(3) long before Altarelli and Cabibo j3| have used them. R.C. Hwa in his 
elaborated work termed them as valon, extended it and showed its application to many physical processes [Q. In 
Rcf.[2] it is suggested that the concept of dressed quark and gluon might be useful in the area of jet physics and 
heavy quark effective theory. Despite the ever presence of CQ no one has calculated its content and partonic structure 
without resorting to hadronic data and the process of deconvolution. The purpose of this paper is threefold: (a) we 
will evaluate the structure of a CQ in the Next-to-Leading Order (NLO) in QCD. (b) we will verify its conformity 
with the structure function data of nucleon and pion for which there are ample data available, (c) In the process, 
however, we will notice that additional refinements are needed to account for the violation of Gottfried Sum Rule 
(GSR) and the effect of binding of CQ's to form a physical hadron. 



II. FORMALISM 



By definition, a CQ is a universal building block for every hadron, that is, its structure is common to all hadrons 
and generated perturbatively. Once its structure is evaluated, in principle it would permit to calculate the structure of 
any hadron. In doing so we will follow the philosophy that in a DIS experiment at high enough Q 2 it is the structure 
of a CQ which is being probed and at sufficiently low value of Q 2 this structure cannot be resolved thus, a CQ behaves 
as valence quark and hadron is viewed as the bound state of its CQ's. Under these criteria partons of DIS experiments 
are components of CQ and at high Q 2 one can write for a U-type CQ its structure as follows: 
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Q 2 ) = \z(fis + Gu ) + l z (G« +G i+ G^+ G») + ... (1) 

where all the functions on the right-hand side are the probability functions for quarks having momentum fraction z 
of a U-type constituent quark at Q 2 . Similar expression can be written for a D-type CQ. Defining the singlet (S) and 
nonsinglet (NS) CQ distribution functions as: 

/ 

G 5 = £(G^+G^) = G / + (2/-l)G M/ ( 2 ) 

i=l 



° NS = E( G ^ ~G±_) = G f - G uf (3) 

i=l 

where Gf is the favored distribution describing the structure function of a quark within a CQ of the same flavor while 
unfavored distribution G u f describes the structure function of any quark of different flavor within the CQ. / is the 
number of active flavors. In terms of singlet and nonsinglet distributions they are as follows: 

G f = ±(G s + (2f-l)G NS ) (4) 

Guf = ±(G s -G NS ) (5) 

Having expressed all the structure functions of CQ's in terms of Gf and G u f we now go to the moment space and 
define the moments of these distributions as: 

M(n,Q 2 )= [ x n - 2 F(x,Q 2 )dx (6) 
Jo 

M 4 (n,g 2 )= f x^GifaQ^dx (7) 
Jo 

the subscription i stands for S, NS. Using charge symmetry, in the following we will refer to CQ distribution only in 
proton. 

In the NLO approximation the dependence of the running coupling constant, a on Q 2 is given by: 

with ft = ±(33 - 2/) and ft = 102 - 2M. The moments of NS and S in the NLO are: 

M NS (n O 2 ) - n + "(Q 2 ) Z a (QQ ( iML. _ A7^ )JV N 1r a.(Q a ) w<SW^ (Q) 
l ' yJ_[ + 4tt [ 2ft 2ft 2 )l[ a s (Q 2 y W 

M s (n,Q 2 ) = {( asiQ2) A[p N - 1 a *(Q 2 o)- a °(Q 2 ) * n * „ , MQ 2 o) _ a J (Q 2 ) ) ^l ) 

a s (Qo) 2ft 47r 47r 47r 

AT N N 

P - 1P - '■(+—)> do) 



2ft 2 + A^-A w - 



where 7^ = — §^7 ( ' ' )Ar and 7^°)^ and are anomalous dimension metrices || and A± denote the eigenvalues 

of one-loop anomalous dimension matrix j(°) N . 
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and p± is given by: 

p£ = ±(7 (0) "-7£)/(A?-A?) (12) 

The quantities d± , the leading order anomalous dimentions, given in Ref.[4], are related to A± in our notations, t is 
the evolution parameter defined as: 

InK 

t = In^r (13) 

The coefficients 7^ : °' )JV and can be found in [||. 

We have taken our initial scale Qq = 0.283 GeV 2 and A = 0.22 GeV. It seems that evolution of parton distributions 
from such a low value of Qo ^ s n °t justified theoretically. The above value of Qo corresponds to a distance of 0.36 fm 
which is roughly equal to or slightly less than the radius of a CQ. It may be objected that such a distances are probably 
too large for a meaningful pure perturbative treatment. We note that F 2 ® \z,Q 2 ) has the property that it becomes 
8{z — 1) as Q 2 is extrapolated to Qq, which is beyond the region of validity. This mathematical boundary condition 
signifies that the internal structure of a CQ cannot be resolved at Qq in the NLO approximation. Consequently, 
when this property is applied to Eq.(19) bellow, the structure function of the nucleon becomes directly related to 
xGcq(x) at those values of Qq, that is, Qq is the leading order effective value at which the hadron can be considered 

as consisting only of three (two) CQ's, for baryons (mesons). In fact our results are only meaningful for Q 2 . > 0.4 
GeV 2 . As it is stated above, the moments of the CQ structure function, \z,Q 2 ) are expressed completely in 
terms of the evolution parameter, t, of Eq. (13). From the theoretical standpoint, both A and Qo depend on the order 
of the moments n. In this work we have assumed that they are independent of n, hence introducing some degrees of 
approximation to the Q 2 evolution of the valence and sea quarks. However, on one hand there are other contributions 
like target-mass effects, which add uncertainties to the theoretical predictions of perturbative QCD, while on the other 
hand since we are dealing with the CQ, there is no experimental data to invalidate an n independent A assumption. 
The moments of valence and sea quarks in a CQ are: 

M^^^M NS {n,Q 2 ) (14) 

M^ = ^(M S -M NS ) (15) 

where M S ' NS are given above. Evaluating Mw^ce and at any Q 2 or t is now straight forward. Using Inverse 

Mcllin Transform techniques, following forms for the valence and sea quark distributions inside a CQ is obtained in 
the NLO: 

zq^i_{z,Q 2 ) = az b {\ - z) c (16) 

CQ 



zq^ (z, Q 2 ) = az\\ - zy[l + V z + £z ' 5 ] (17) 

The parameters a, 6, c , a, etc. are functions of Q 2 through the evolution parameter t. The same form as in Eq.(17) is 
obtained for Gluon distribution in a CQ but only with different parameters. Functional form of them is a polynomial 
of order three in t and are given in the appendix. We notice that the following sum rule reflecting the fact that each 
CQ contains only one valence quark is satisfied for all values of Q 2 : 

J q^_{z,Q 2 )dz = l. (18) 

Substituting these results in Eq.(l) completes the evaluation of a constituent quark structure function in NLO. In 
Figure (1) various parton distributions inside a CQ is plotted. 
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III. HADRONIC STRUCTURE 



In previous section we calculated the NLO structure of a CQ. In this section we will use the convolution theorem 
to calculate the structure function of proton, F^x, Q 2 ), and that of a pion. Let us denote the structure function of 
a CQ by F 2 (z, Q 2 ) and the probability of finding a CQ carrying momentum fraction y of the hadron by Gcq (y). 

h 

The corresponding structure function of the hadron, using the convolution theorem is as follows: 

F 2 h (x,Q 2 )=J2 / 1 -Gc £ (y)^(^,Q 2 ) (19) 

CQ J * V " V 

where summation runs over the number of CQ's in a particular hadron. Also notice that Gco_{y) is independent of the 
nature of the probe and its Q 2 value. Gcq (y) in effect, describes the wave function of hadron in CQ representation 

h 

containing all the complications due to confinement. From the theoretical point of view this function cannot be 
evaluated accurately. To facilitate phenomenological analysis, following Ref.[4] we assume a simple form for the 
exclusive CQ distribution in proton and pion as follows: 

Guud/ p {Vi,V2,V3) = l(yiy2) m y%S( yi +y 2 + y 3 - 1) (20) 

God/,- (2/1,2/2) - qy?y%6( yi +J/2-1) (21) 
Integrating over unwanted momenta, we can arrive at inclusive distribution of individual CQ: 

G °^ = B{b+1 ] 2a+2 / (i-y) 2a+1 ( 23 ) 
G ^ = wrbrT) r{1 - yr (24) 

similarly expression for Gd/tt- with the interchange of fi <-> v. In the above equations B(i,j) is Euler Beta function 
and its arguments are fixed using the sum rule: 



L 



i 

Goa(y)dy = l (25) 



where CQ = U,D,U and h = p, tt~ . Numerical values are: fi = 0.01, v = 0.06, a = 0.65 and b = 0.35. In Figure (2) 
the CQ distributions in proton and tt~ are shown. We stress that CQ distributions in hadrons are independent of Q 2 
and the nature of probe being used. Now it is possible to determine various parton distributions in a hadron. For 
proton we can write: 

f dy f ^ d y 

<lvai./p(x,Q 2 ) = 2/ — G u/p (y)q vaL/u (x,Q 2 ) + / — G D / p (y)q vaL / D (x 1 Q 2 ) = u vaL / p (x, Q 2 ) + d va i./ p (x, Q 2 ) (26) 

Jx y J x y 

q S ea/ P (x,Q 2 ) = 2 / —G u/p (y)q sea/u (x,Q 2 )+ —G D/p (y)q sea/D {x,Q 2 ) (27) 

J x y J x y 

The above equation represents the contribution of constituent quarks to the nucleon sea. Comparing with data on 
proton structure functions shows that the results fall short of representing the experimental data by a 3-5 percent. 
This is due in our opinion, to the fact that CQ is not free in a hadron but they interact with each other in forming the 
bound states. That means, there are some soft gluons in the nucleon besides the CQ's. In the process of formation of 
bound state, a CQ emits gluon which in turn decays into q — q pairs which gives a residual component to the partons 
in a hadron. In our picture there is no room in the CQ structure for breaking the SU (2) symmetry of the sea but after 
creation of q — q pair from the emitted gluon , these quarks can recombine with CQ to fluctuate into meson-nucleon 
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state which breaks the symmetry of the nucleon see. In what follows we will compute this component, its contribution 
to the nucleon structure function and violation of Gottfried sum rule following the preswcription used in ||. In 
order to distinguish these partons from those confined inside the CQ, we will term them as inherent partons. Although 
this component is intimately related to the bound state problem, and hence it has a non-perturbative origin, for not 
so small values of Q 2 we will calculate it perturbatively for the process of CQ — > CQ + gluon — > q — q at an initial 
value of Q 2 = 0.65GeV 2 where a s is still small enough. The corresponding splitting functions are as follows: 



2 

For the joint probability distribution of the process at hand, we get 



P qg (z) = l:(z 2 + (l + z) 2 ) (29) 



qinhA ,,()^ = qinh (x,Q 2 )=N 7 ^ t—Pisi-) t -P 9q (-)G CQ (z) (■■!()! 

v^"/ Jx y y Jy z z 



The splitting functions and the qi n h(x,Q 2 ) = q~inh{x,Q 2 ), above are that of the leading order rather than NLO. We 
do not expect it should make much of a difference since, its contribution to the whole structure function is only a 
few percent as can be seen in Figure (3). In the above equation N is a factor depending on Q 2 and Gcq is the 
constituent quark distribution in the proton given previously. The same process, however, also can be a source of 
SU(2) symmetry breaking of nucleon sea and resulting in u sea ^ d sea , and hence the violation of Gottfried sum rule 
(GSR). There are several explanations forthis observation such as flavor asymmetry of the nucleon sea || ]l0[ ], isospin 
symmetry breaking between proton and neutron, Pauli blocking, etc. One of these explanations fits well within our 
model. It was proposed by Eichten, Inchliffe and Quigg |Tl| that valence quark fluctuates into quark and a pion. 
In other words a nucleon can fluctuate into a meson-nucleon state. This idea is appealing in our model and can be 
calculated rather easily. In our model after a pair of inherent q — q created a u can couple to a D-type CQ to form 
an intermediate n~ — Du while the u quark combines with the other two U-type CQ's to form a A ++ . This is the 
lowest uu fluctuation. Similarly a dd can fluctuate into the 7r + n state. Since A ++ state is more massive than n state, 
then the probability of dd fluctuation will dominate over uu fluctuation which naturally leads to an excess of dd pairs 
over uu in the proton sea. This process is depicted in Figure (4). Probability of formation of a meson-barion state 
can be written as in Ref . [8] : 

Pmb(x)= f —f — F(y,z)R(y,z;x) (31) 

Jo y Jo z 

where F{y 1 z) is the joint probability of finding a CQ with momentum fraction y and an inherent quark or anti-quark 
of momentum fraction z in the proton. R{y, z; x) is the probability of recombining a CQ of momentum y with an 
inherent quark of momentum z to form a meson of momentum fraction x in the proton. The evaluation of both of 
these probability functions are discussed in [l2| for a more general case and an earlier, but pioneering, version also 
proposed in [pUf . In the present case these functions are much simpler. Guided by works done in Ref. [8,12,13] we 
can write : 

F(y,z) = n y GcQ(y)zq mh .(z)(l-y-z) s (32) 

P 

R(y,z;x) = py a z b 5(y + z-l) (33) 

Here we take a = b = 1 reflecting that two CQ's in meson almost equally share its momentum. The exponent S is 
fixed for the n and A ++ states using the data from E866 experiment and the mass ratio of A to n. They turn out 
to be approximately 18 and 13 respectively. f2 and p are the normalization constants also fixed by data. It is now 
possible to evaluate um and du quarks associated with the formation of meson states: 

d M (x,Q 2 ) - f 1 ^[PUy) + 7^a++(2/)K(-,Q 2 ) (34) 
Jx y 6 y 

1 f 1 dy x 

um(x,Q 2 ) = - — f7rA++(2/)w7r(-,Q 2 ) (35) 

2 J x y y 
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where and are the valence quark probability densities in the pion at scale Qq. The coefficients | and g are due 
to Isospin consideration. Using Eqs. (16, 17, 24) we can calculate various parton distributions in a the pion. Those 
pertinent to Eqs. (34, 35) are: 

i2n _ f 1 ^ — _, x ^2\ d V 



<ai.M)= / Go^y^^oi-^j (36) 



< Q ,(x,Q 2 )= / G Dh -{y)d vaL , D {-,Q 2 )^- (37) 

Jx y v 

Uval./U = Uval./U (38) 

There are some data on the valence structure function of 7r _ [ jl4| . Defining valence structure function of 7r _ as: 

^roii. = X ^val. — xd val. (39) 

we present the results of our calculation for F* al in Figure (5) along with the experimental data at Q 2 around 6 Gev 2 . 
Returning to the d excess over u in proton we can write all the contributions as: 

dwo„ = fr + fo-+cq (40) 

U u M + U inh . +C Q 

NuSea collaboration at FermiLab E866 experiment has published their results |H| for integral of d — u and i at 
Q = 7.35 GeV. With the procedure described, we have calculated these values at the same Q and for the range of x 
as experimentally measured: x — [0.02,0.35]. The results of the model is: 

0.345 

dx(d - u) = 0.085 (41) 

.02 

to be compared with the experimental value of 0.068 ± 0.0106. We get for the entire range in x: J^dx(d~u) = 0.103 
while the experimentally extrapolated value is 0.1 ± 0.018 which are in excellent agreement with our calculations. 
This gives a value for the Gottfried sum rule of: 

P 1 dr 1 If 1 

S G = / [F*(x) - F?(x)]— = - - - / dx[d(x) - u(x)} = 0.264. (42) 

JO X O 6 Jo 

at Q = 7.35 GeV. The NMC result @ is S G = 0.235 ±0.026 which is at much lower value of Q 2 = 4 GeV 2 . In Figure 
(6), d(x) — u{x) and 4 in proton are shown as a function of x at Q — 7.35 GeV along with the measured results. 
We are now in a position to present the results for the proton structure function, F%. In Eqs. (16, 17) we presented the 
form of parton distributions in a CQ. Using those relations with the numerical values given in appendix then from 
Eq.(l) the structure function of CQ is obtained. In Eqs. (22, 23) the shape of CQ distributions in proton is given. 
With the help of Eq. (19) now all the ingredients are in place to calculate F%. In Figure (7) the results are shown 
at many values of Q 2 . As it is evident they fall a few percent short of representing the data. However, as mentioned 
earlier, there is an additional contribution from the inherent partons to F% which is calculated in Eq. (30). Adding 
this component represents the data rather well and can be seen from Figure (7). The data points are from Jl7|| . For 
the purpose of comparing our results with other calculations, we have also included in Figure (7), the GRV's NLO 
results Q as well as the prediction of CTEQ4M |@. Notice that we have taken the number of active flavors to 
be three for Q 2 < 5GeV 2 and four flavors elsewhere. In Figure (8) the gluon distribution predicted by the model is 
presented along with those from Ref. [18, 19]. 

IV. SUMMARY AND CONCLUSION 

In this paper we have used the notion of constituent quark as a well defined entity being common to all hadrons. Its 
structure can be calculated in QCD perturbatively to all orders. A CQ receives its own structure by dressing a valence 
quark with gluon and q — q pairs in QCD. We have calculated its structure function in the Next-to-Leading order 
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for the first time. Considering a hadron as the bound states of these CQs we have used the convolution theorem to 
extract the hadronic structure functions for proton and pion. Besides the CQ structure contribution to the hadrons, 
there is also a nonperturbative component while contributing only a few percent to the overall structure of hadrons 
becomes crucial in explaining the violation of Gottfried sum rule and the excess of d over u in the nucleon sea. A 
mechanism is devised for this purpose and necessary calculations are outlined. We have presented the results and 
compared them with all available relevant data and with the work of others. We found that our results are in good 
agreement with the data. 



V. APPENDIX 



In this appendix we will give the functional form of parameters of Eqs. (16, 17) in terms of the evolution parameter, 
*. This will completely determines partonic structure of CQ and their evolution. The results are valid for three and 
four flavors, although the flavor number is not explicitly present but they have entered in through the calculation of 
moments. As we explained in the text, we have taken the number of flavors to be three for Q 2 < hGeV 2 and four for 
higher Q 2 values. 

I) Valence quark in CQ (Eq. 16): 

a = -0.1512 + 1.785* - 1.145* 2 + 0.2168* 3 
b = 1.460 - 1.137* + 0.471* 2 - 0.089* 3 
c = -1.031 + 1.037* - 0.023* 2 + 0.0075* 3 

II) Sea quark in CQ (Eq. 17): 

a = 0.070 - 0.213* + 0.247* 2 - 0.080* 3 
(3 = 0.336 - 1.703* + 1.495* 2 - 0.455* 3 
7 = -20.526 + 57.495* - 46.892* 2 + 12.057* 3 
?] = 3.187 - 9.141* + 10.000* 2 - 3.306* 3 
£ = -7.914 + 19.177* - 18.023* 2 + 5.279* 3 
N = 1.023 + 0.124* - 2.306* 2 + 1.965* 3 

III) Gluon in CQ (Eq. 17) 

a = 0.826 - 1.643*+ 1.856* 2 - 0.564* 3 
f3 = 0.328 - 1.363* + 0.950* 2 - 0.242* 3 
7 = -0.482 + 1.528* - 0.223* 2 - 0.023* 3 
j] = 0.480 - 3.386* + 4.616* 2 - 1.441* 3 
£ = -2.375 + 6.873* - 7.458* 2 + 2.161* 3 
N = 2.247 - 6.903* + 6.879* 2 - 1.876* 3 
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VI. FIGURE CAPTION 



Figure- 1. Moments of partons in a CQ at Q 2 = 20 GeV 2 as a function of z. 
Figure-2. Parton distributions in proton and tt~ at Q 2 — 20 GeV 2 as a function of x. 

Figure-3. Contribution of inherent component to a sea parton distribution. The dashed-dotted line is that of CQ and 
the solid line represents the sum of the two components. 

Figurc-4. Processes responsible to SU(2) symmetry breaking in the nucleon sea and violation of Gottfried sum rule. 
Figure-5. Pion valence structure function as a functuon of x at Q 2 = 5.5 (Gev/c) 2 . Solid curve is the result of the 
model calculations and the data points are from Ref.[12]. 

Figurc-6. The ratio 4 and the difference d — u as a function of X. The solid line in the model calculation and the 
dotted line is the prediction of CTEQ4M. Data are from Ref. [13]. 

Figure-7. Proton structure function F$ as a function of x calculated using the model and compared with the data 
from Ref. [15] for different Q 2 values. The thin line is the prediction of GRV Ref. [16] and the dashed line is that of 
CTEQ4M Ref. [17]. 

Figure-8. The gluon distribution in proton as a function of x at Q 2 = 20 GeV 2 . We have also shown the prediction 
of GRV (dashed-dotted line) and CTEQ4M (dashed line). The data points are from HI collaboration. 
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Fig. 7 (continued) 
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